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Chapter 1 



GRAVITATIONAL WAVES 

1.1 Introduction 

Almost every student of science is familiar with the phenomenon of waves, such 
as water waves ( ripples rolling across the ocean) , sound waves (vibration in the 
air), etc. The above mentioned kinds of waves are easily observable. Electro- 
magnetic waves are comparatively difficult to understand but not as difficult 
as gravitational waves. The existence of gravitational waves was disputed for 
a long time but is generally accepted now. 

What are gravitational waves? How do they propagate? and what is their 
energy content? These questions are addressed in the first two chapters. In the 
third chapter the pseudo-Newtonian formalism and its extension is reviewed 
in general and the formula for the momentum imparted to test particles in 
arbitrary spacetimes is reviewed in particular. In chapter four the analysis of 
a paper claiming to determine the spin for gravitational waves is given, and 
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compared with the spin given by a geodesic analysis. It is demonstrated that 
the other claim is inconsistent. Finally in chapter five a summary of the work 
is given with the conclusion. 



1.2 



Idealization 



The only way to come to grips with so complicated a subject as General Rel- 
ativity is by idealization . Study one idealization after another. Build a cata- 
logue of idealization, of their properties and of techniques for analyzing them 
and then arrive at a conclusion. 

Let's now see how can we idealize gravitational waves. As one idealizes 
"water waves" as small ripples of geometry rolling across the ocean so one 
gives the name "gravitational waves" to small ripples rolling across the space 
time. Both of the waves are idealizations. One cannot, with infinite accuracy, 
delineate at any moment which drops of water are in the waves and which 
are in the underlying ocean. Similarly one cannot tell precisely which parts of 
spacetime are in the ripples and which are in the cosmological backgrounds. 
One can almost do so; otherwise one would not speak of "waves" . Look at the 
ocean, the seascape is dominated by the waves. Changes occur at the surface 
of the ocean, which propagates obeying the following [1] wave equation: 



Similarly gravitational waves are perturbations of spacetime. As I have already 




1 92 92 



) (height of the surface) = 0. 




g"^ dt^ dy'^ dx^ 
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stated that gravitational waves are ripples of geometry. So let me support my 
statement. Take a massive body and disturb it violently, the near field adjusts 
rapidly, but the far field must wait for the signal that the mass has moved to 
propagate to it at a finite speed "c" . Thus there is a travelling kink which falls 
off in strength with distance. Hence the gravitational waves are small ripples 
rolling across the spacetime. 

Now get more sophisticated. Notice from a space ship the large-scale cur- 
vature of the ocean's surface-curvature because the Earth is round. Curvature 
because the Earth, Sun and Moon pull the water. As waves propagate long 
distances, this curvature bends their fronts and changes slightly their simple 
wave equation. Spacetime is similar. Propagating through the universe, ac- 
cording to Einstein's theory must be a complex pattern of small-scale ripples 
in the spacetime curvature, ripples produced by binary stars, by gravitational 
collapse, by explosion in galactic nuclei etc. 

1.3 Linear approximation for the investiga- 
tion of gravitational waves 

We are discussing gravitational waves in the frame work of General Relativity, a 
nonlinear field theory of gravity. Though many of the interesting consequences 
of General Relativity comes from its non linearity it is worthwhile to study its 
linear approximation. Linearization actually leads one to the gravitational 
waves. 
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What makes General Relativity nonlinear? As the field in General Rela- 
tivity is the metric tensor, its appearance in the field equations non-linearly 
gives rise to the non-linearity of General Relativity. We cannot change the 
way the metric tensor enters into the curvature but we can write the curved 
spacetime metric as the fiat spacetime metric tensor, 77^,^, and an additional 
term, hf^^- We then require that hf^u and its derivative occur only once in the 
field equations and higher powers be neglected [6]. 



(1.2) 



Let 



9' 




(1.3) 



where rj^'' is the inverse of 77^1^. Then by definition. 



9'%.-{r'+fn{vp.+K.), 



(1.4) 



Which implies that 



(1.5) 



Cancelling 5^ on both sides and multiplying through by rf^ we get 



(1.6) 
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The last term is clearly quadratic in the difference between the curved and flat 
spacetime metric tensor. Thus to first order, 

/M-^_^MP^-/i^^ + 0(/i2). (1.7) 

Using the flat spacetime metric tensor to raise and lower indices we can write 

= -h"" + 0{h^). (1.8) 
Using equation (1.3), equation (1.8) becomes 



Using this linearization and for the moment taking Cartesian coordinates, 
so that there are no derivatives of 77'*'^,the Christoffel symbols linearize to: 



/i V 

Clearly terms quadratic in the Christoffel symbols become quadratic in h and 
can be neglected compared with linear terms. Thus the linearized Ricci tensor 
is 




~ 'i^^l'^ {hixTT,!/ -\- huT^^n — hnu^Tj)\^p (-'-■-'-2) 

2^^^ (^/U7r,i^p ~l~ ^i/7r,/ip h^ii^pj^ hp^^nn^ (1.14) 

A choice of coordinates can be made to have the first two and the last term 
in the brackets disappear. To see this first we note that we can rewrite the 
Ricci tensor as: 



-\r]'^h„, (1.15) 

where h — hf, then we can break the first bracket in to two terms. Now consider 
an infinitesimal transformation 

xi" ^x'^^x^ + ^^ix"), (1.16) 

so that the terms quadratic in /i or ^ can be neglected. Since this is only a 
coordinate transformation, it must leave the metric invariant and hence 



ds'^ = g^,r,{xP)dx>'dx'' ^g'^^{x'P)dx'i'dx"' 

= g'^^{x"^){dx^-ri%dx'^){dx^-ri:pdxf'). (1.17) 

Using the linearization procedure it is easy to see that 

^MZ. ~ V - ^M,^ - (1-18) 
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Thus we have 



-C. + C*^^- (1-19) 

Differentiating the above equation relative to x'^ it is easy to see that the 
last two terms simply cancel and we have 



(K - = - ^^'^P.P - ^''^M,.P- (1-20) 

Taking the harmonic gauge condition 



we can make 



ih', - ^M^,, = r,,p = 0. (1.22) 

It should be pointed out here that the use of Cartesian coordinates was not 
crucial, but merely for convenience. With any other coordinates we could have 
to introduce the corresponding flat spacetime Christoffel symbols. This has 
been avoided so as not cause confusion of notation. We now drop the primes 
and so obtain 

R^.u ~ -^r?^"V,P- = -^DV- (1-23) 

We can contract Equation (11.221) to obtain the Ricci scalar {R = —^nh) 
and hence combine them to obtain the Einstein tensor. Thus the Einstein field 
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equation becomes 

□0^, = -2kT^,. (1.24) 

Here T^,^ is the stress energy tensor and k is the proportionahty constant. In 
regions where T^i^ = 0, (f)^^ then satisfies the wave equation, Thus (f)fj_v represent 
gravitational waves. 

1.4 Plane wave solution in the linearized the- 
ory 

The simplest of all solutions to the linearized equation /i^j, q, = is the monochro- 
matic plane wave solution [1] 

V = ^[Vexp(^A;„a;")]. (1.25) 

Here 3? means that one must take the real part in the bracket. While A^j^^, is 
the amplitude and is the wave vector, satisfying 

kak"" = (k a null vector), (1.26) 
■^nak'^ — (A orthogonal to k). 
This solution describes a wave with the frequency 

^^k'^{ki + ki + kiy/', 

which propagate with the speed of light in the direction {l/k^){kx, ky, kz). At 
first sight the amplitude Aj^^, appears to have six independent components 
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(ten, less the four orthogonality constraints ^^a^" — 0). But this cannot be 
right. The gravitational field has two dynamic degrees of freedom (not six). 
The plane wave vector 



with four arbitrary constants C^, generates a gauge transformation that can 
change arbitrarily four of the six independent components of ^4^,^. One gets rid 
of this arbitrariness by choosing a specific gauge. 

1.5 The transverse traceless (TT) gauge 

Consider the four velocity u through all the space time and impose the condi- 
tion [1] 



These are only three constraints on Aj^^. Why not four? Because A is orthog- 
onal to k i.e. k^^A^j^i^u") = 0, as already mentioned. As a fourth constraint 
use a gauge transformation to set A/j^ — 0. So there are now four constraints in 
all, Aj^aU" = A^k'^ = A^ = 0, on the ten components of amplitude. Thus the 
two remaining free components of A^^, represents the two degrees of freedom 
in the plane gravitational wave. 

It is useful to restate the four constraints, A^u^ = A^ak^ = A^ = in the 
Lorentz frame where vP — l,u^ — Q and in a frame where /c" does not appear 



^'^^ -iC^exp(iA;«x°), 



(1.27) 




(1.28) 
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explicitly: 



hi_io = , i.e only the spatial componemts hjk are non-zero (1-29) 
hkj,j = i.e the spatial components are divergence free (1.30) 
hkk = = i.e the spatial components are trace free. (1-31) 

The gauge conditions are all linear in /i^,^ so an arbitrary wave will also satisfy 
the above gauge conditions. In this gauge only the hjk are non-zero. So we 
need only to impose the six wave equations 

Dhjk = hjk^aa = (1.32) 

Now the definition. Any symmetric tensor satisfying the constraints (11.291) . 
(I1.30p and (11.311) is called a transverse traceless tensor. Why it is called trans- 
verse traceless? Because it is purely spatial {h^o = 0) and, if thought of as 
a wave is transverse to the direction of propagation {hijj = hijkj = 0) and 
traceless because hkk = 0. 

The special gauge in which /i^,^ reduces to its transverse traceless part is 
called the transverse traceless gauge . The conditions (11.29p . ( 11.30P and (]1.3ip 
defining this gauge can be summarized as: 

V = (1.33) 

Only pure waves can be reduced to TT gauge. In the TT-gauge the space 
components 

RjOko = RojOk = —RjOOk = —RojkO- (1-34) 
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of the Riemannian curvature tensor have an especially simple form 




(1.35) 



As the curvature tensor is gauge invariant therefore hf^^, cannot be reduced to 
still fewer components than it has in the TT-gauge. 

1.6 Comparison with electromagnetic waves 

A simple system consisting of two charges of equal magnitude but of opposite 
signs, each situated at a distance r/2 from the origin O (see Fig.1.1), which is 
taken to lie on the line connecting the charges, is the simplest example of an 
electric dipole. The oscillation in the dipole generates electromagnetic waves. 

Now we define the electric dipole moment of the pair of equal charges as 
the product of charge q and the separation r. 



where e is the unit vector from negative to positive charge. Consider a system 
of charges Qa and let Va be their radius vectors then 



is called the dipole moment of the system of charges. 

Consider masses in an isolated, nearly Newtonian system, moving about 
each other. They emit radiations. For an order of magnitude estimate, one 



d = qre, 



(1.36) 





a 
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Figure 1.1: A system consisting of four masses attached by springs, forming 
two dipoles. The oscillation in the system generate gravitational waves. 
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can apply the familiar radiation formula of electromagnetic theory, with the 
replacement —>■ —rn? which converts the static Coulomb's law into Newton's 
law of attraction. Although it introduces a moderate error in numerical factor 
and changes angular distribution, but it gives an estimate of the total power 
radiated. In electromagnetic theory, electric dipole radiation dominates, with 
power output or luminosity "L" given by[l], 

^elecdip = g*?^*^^' (1.38) 

where a is the acceleration in the dipole. Then using equation fll.36p we can 
write 

2 -2 

Le.ec..p =-d . (1.39) 

The gravitational analogue of the electric dipole moment is the mass dipole 
moment, 

d = '^rriAXA, (1.40) 

A 

its first time-rate of change is the total momentum of the system, 

d=J2^AXA=P, (1.41) 

where P is the momentum of the system. The time-rate of change of the mass 
dipole moment has to vanish because of the law of conservation of momen- 
tum, d = P = therefore there can be no mass dipole radiation in gravitation 
physics. 
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The next strongest type of electromagnetic radiation are magnetic dipole 
and electric quadrupole. Magnetic dipole radiation is generated by the sec- 
ond time derivative of the magnetic moment, /j^ . Here again the gravitational 
analogue is a constant of motion of the angular momentum i.e., 



fj, — y {position of A) x {current due to A) 

= '^ta X {mVA) = J = constant (1-42) 

Which shows that /i= 0, so it can not radiate. Thus, there can be no gravita- 
tional dipole radiation of any sort. Physical example is a system of two masses 
attached by a spring, acts as an oscillating dipole. 

1.6.1 Comparison with plane electromagnetic waves 

Consider the metric [1] (with signature — , — , — )) 

ds^^L\u){dx^ + dy^)-dudv. (1.43) 

where u — t — z andv — t + z which is always flat. It satisfies the vacuum 
Einstein equations R^i, = which implies that L" = (see Appendix 1 for 
the proof). Now if L" — the spacetime is static therefore cannot represent 
gravitational waves. In this metric the electromagnetic potential 

A = A^dx^ = A{u)dx (1.44) 

satisfies the Maxwell equation for arbitrary A{u). It represents an electromag- 
netic plane wave analogous to the gravitational plane wave. The only non-zero 
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components of this wave are 



Fu^ = A i.e Ft. = -F,. = A'. (1.45) 

So the wave propagates in the z-direction, the magnetic vector oscillates in 
the y-direction and electric vector in the x-direction. The only non-zero vector 
of the stress energy tensor is 

T„„ = (47rL2)-i(A')^ (1.46) 

It can easily be verified that the Maxwell equation are satisfied by fll.44p in 
dUS]). 

To make the metric acceptable, we need to impose the Einstein equations 
Rfiu — \gnvR = ^^^41^5 ^ is generally ^^,but here we use gravitational units i.e. 
G = c = 1 so it becomes Svr. As all the nonvanishing Ricci tensor components 
are the same and they read , so the Einstein field equations becomes 
L" + [AttT^^L = 0. This has exactly the form of the equation L" + {(3'YL = 0, 
(which will be discussed in detail in chapter 2) for gravitational plane wave. 
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Chapter 2 

REALITY OF 
GRAVITATIONAL WAVES 

In the linearization procedure the stress energy tensor was taken to be zero, 
which creates a conceptual problem, namely the question of reality of gravita- 
tional waves. This will be discussed in detail in the first section of this chapter. 
For this purpose some exact gravitational wave solutions of the Einstein field 
equations are presented. Returning to the purpose of this chapter the energy 
contents of the waves are discussed. Finally as an additional but necessary, 
topic some sources of gravitational waves are discussed with in the scope of 
this work. 
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2.1 The conceptual problem 

The exact gravitational waves being by definition, solutions of the vacuum field 
equations, have a zero stress energy tensor. This creates a conceptual problem. 
How can these solutions of the Einstein field equations represent waves if they 
carry no energy? Essentially, this problem arises because energy is not a well 
defined concept in General Relativity. For energy to be well defined in General 
Relativity the metric must have a time like isometry (Killing vector) so as to 
allow time translational invariance. This will not generally be true. Infact, 
spacetimes for which it is true are static whereas gravitational wave solutions 
must be non-static. Thus energy is not well defined for spacetimes containing 
gravitational waves. The question then is, how can we tell that these solutions 
really do behave as we would expect of the waves? 

One way to answer this question: the very process of linearization provides 
the energy. The point is that if Dcf)^,^ — exactly then R^jy to order h 
but 7^ to higher orders in h. Thus T^i, ^ only to order h but T^i, ^ 
generally. Conversely if T^y = and hence = exactly then D^^y ~ 
only to order h. Thus we can expand Rf^u in powers of h, retaining linear terms 
on the left side of the equation and transposing all higher powers to the right 
side. These higher order terms become an effective stress energy tensor and 
the linearized equations give the gravitational waves. 
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2.2 Some exact solutions of gravitational waves 

So far we have obtained the wave equations for gravity by hnearizing the 
Einstein field equations. In principle the solutions so obtained could be exact 
solutions of the vacuum Einstein field equations, of course there could be trivial 
static solutions which effectively satisfy the Laplace equation. But they do not 
represent moving waves so we are interested in the solutions which are non- 
static. 

The first solution to be discovered was for cylindrical gravitational waves, 
by Einstein and Rosen in (1937) [2]. So first consider this: 

2.2.1 Cylindrical gravitational wave solution 

A cylindrically symmetric metric depending on two arbitrary functions 7 and 
■0 of the time t and cylindrical radial coordinate p, is 



(2.1) 



The metric tensor is 



Qab 



,2(7-V) 










-p^e-^V- 




(2.2) 



Its inverse is 
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g-2(7-V) 

-e-2(7-V') 








non-zero Christoffel symbols are: 



V 








-e^^ 





hi 




hi 


\ ^ 1 


hi 




[ooj 




[-J 











7' --0 





hi 




hi 


\ ^ 1 


hi 




[o.J 




i-J 




[ooJ 




[-1 




2 2 




3 3 



, ^^.g2(2V,-7). 
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2 2 
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[3] 


1-^ 1 
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hi 








[30J 




i-1 
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H 



(inv^),o = 2(r-v^-); 



(2.9) 



fx 

fXl 



^ = (lnV^),i = 2(7'-V'') + l/p. 



(2.10) 



The ""' refers to differentiation with respect to p and the dot differentia- 
tion with respect t. The non vanishing components of the Ricci tensor are 
Rm-i Ru, R22 and it^sa.Only four components are needed for the purpose. Us- 
ing the Christoffel symbols hsted above we obtain the following Ricci tensor 
components, 



■00 



A* I I A* I I A* I I ^ 

> -(lnv^),oo + (lnv^),^ 

0// |00||Oz/||0/i 

> +< > - (liiV9),oo+ (lnv^),o 










^ ] 






2 




[osj 



(2. 
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After simplification we get 



i?oo = r-r + i'-r-'^{r-r) + {i-iJ-?-'^{i-i^'f 
-a - ^-f - 2(^f + [1 + 2(y - ^o](7' - ^0 

= -(7" - r) + (7" - r) + -(V - V'O - 2(V'-)'- (2.13) 

p 

Similarly other non-vanishing components are, 



Rii = (7 ■ - - (7" - + -(7 + - 2{rf. (2.14) 



p 

(2.15) 

i?oi = -1 - 2V'-^'. 
P 

(2.16) 
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So the Einstein field equations becomes 



(7- - r) - (7" - - -(7 - + 2(^-)' = 0. (2.17) 

P 

(7- - - (7" - + ^(7' + - = 0. (2.18) 

P 

(2.19) 

-7--2^-?/'' = 0. 
P 

(2.20) 

Equation (12.191) is the usual cylindrical form of the wave . As this is a second 
order linear differential equation, the solution has two arbitrary constants, one 
corresponding to the ingoing cylindrical wave and the other corresponding 
to the outgoing. Retaining only the outgoing waves with amplitude A and 
frequency uj we have 

?/;(t,p) = A[Joix)cos{ujt) + No{x)sm{ujt)]. (2.21) 

where x = up, Jo{x) and No{x) are the zero order Bessel and Neuman func- 
tions respectively. Add equation fl2.17l) and fl2.18p to obtain 
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7 = P{f^ + (2.22) 

Equations (12.201) and fl2.22p gives the space and time derivative of 7(t, p) in 
terms of functions that are now known through equation fl2.2ip . The only thing 
required is the integration with respect to each variable. The time integration 
is relatively easy and the space integration, though tedious, is in principle 
easy (using the standard formulae for the integrals of the Bessel and Neumann 
function). The resulting solution of 7 is then 

7(t,p) = \A^x{Mx)J'^{x) + No{x)N'^{x)+x[J'^{xf + N'^{xf] 

+ [Jq{x)J'q{x) - Nq{x)N'q{x)] cos{2ujt) ) ■ (2-23) 

+ [Jo{x)N'^{x) + J'q{x)Nq{x)] sin(2cjt)} - ^A^uJi 
Where prime now refers to differentiation with respect to x not p. Hence the 
metric 

ds^ = e2(^-^)(dt2 _ dp') _ e- - e^'^dz'. (2.24) 
represents cylindrical gravitational wave with the above definition of 7 and ip . 

2.2.2 Plane gravitational wave solution 

The solution which we are going to present here was discovered by Bondi 
and Robinson in 1957 [3]. We take a line element which incorporates the 
symmetries of a plane and represents a wave going in the x-direction. 
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ds^ = e^^^^\de - dx") - u\e^^^^^dy^ - e"^^^^^dz^). (2.25) 

Here all coefficients in the metric are functions of (t— x) which are represented 
by u. The (x, z) are not the usual Cartesian coordinates but rectangular 
coordinates in a curved space-time. Thus the metric tensor of equation fl2.25p 
is 



9ah 



g2n(«) g Q 

-e'^^i^^ 
-u2e2/^(") 




Its inverse can be written as 
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(2.26) 
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Here it should be noted for convenience that 

Vtfl = Q! = — and (3^ = [3' = — Here " / " represents derivative with 
respect to u. 
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The non zero Ricci tensor components are Rqo, Rqi andi?ii. 



^00 = <\ } -(In v^), 00+ (In v/^),^ 



/i 



f ^ 1 


hi 


\ . 1 
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I ^ \ 


[ooj 








[ /i J 



(2.32) 

For the given metric (12.251) we have 



(In v^^) =2\nu+2Q. 



(2.33) 
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Inv1^)' = 2(fi' + i (2.34) 
u 



(In v^)" = 2{n" - 1) (2.35) 
Now using the Christoffel symbols listed in equation fl2.3ip we get. 
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Roo = An'/u-2p'^. (2.36) 

Here the vacuum Einstein equation reduces to 



u 



It is an exact plane gravitational wave solution. 



Roi = -4— + 2/5'2 = 0. 
u 
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(2.37) 



(2.38) 



R^^ =4 2/5' 



1/2 



0. 



u 



(2.39) 



Similarly equations (12.381) and (I2.39P reduces to same value as equation (12.371) , 
where other components of Ricci tensor are zero. Hence the only non-trivial 
exact solution of the plane gravitational wave is, 



2.3 Interaction of a particle with plane grav- 
itational waves 

Here it will be shown that gravitational waves carry energy following the 
method of Weber and Wheeler [3], considering the plane gravitational wave 
(I2.25P , which is discussed in the previous section. 

As already mentioned, here we will show that plane gravitational waves 
carry energy. It can be shown by analyzing the motion of a particle which 
is initially at rest and interacts with the plane gravitational wave. Write the 
geodesic equation. 




(2.40) 




(2.41) 
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Using the Christoffel symbols, listed in (12.311) . and the usual summation con- 
vention over the repeated indices, equation (12.411) becomes: 

+ n'i^y + 2fi'$^ - n'i^y + uiu(3' + i)e2{/3-^)(^)2 
as-^ as as as as ay 

as 

(2.42) 

Here in the above equation (12.421) ^ is unknown. It can easily be found as 
follows. Dividing both sides of (12.251) by ds"^ and simplifying, we get 



^ = [(1 + (-2fi + - + ... + (^)2]Vl (2.43) 
as o as 

Using the binomial expansion, we get the following after further simplification. 



^ = i_n + n'-ln' + ... + A' + ... (2.44) 

as 3 as 

Consider as a first order quantity and /3 a second order quantity. Let 

= a;(o)+a;(i)+X(2) + ... (2.45) 

Using these considerations and equation (12.441) . equation (l2.42l) gives the fol- 
lowing approximation equations. 

d'^X/n) ^, dt dX/n) 



+ 2Q'^^-^ - 2f]'^^ = 0. (2.47) 

ds'^ ds ds ds ds 
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^'^(2) I on/(^^ dx(o)^dx^2) n^Jj^^^'^f -Q (2 48) 

ds'^ ds ds ds ds 

Equations (12 .46^ . (12.471) and (12.481) are zero order, first order and second 
order approximation equations respectively. With the previously developed 
tools we can integrate the approximation equations. Here it should also be 
noted that as is a first order quantity, its derivative will also be a first order 
quantity. The result is: 

i) The zero order approximation 



dx 



(0) 



ds 

ii) The first order approximation 



-2n'{l-n)x(o)+Ci, (2.49) 



^^(1) - -2c2n'x^i) + cs; (2.50) 



ds 

in) The second order approximation 



dX(2) 



-2c4^'x^2) + 2c^Vl's + ce. (2.51) 



ds 

Thus equations (12.491) . (I2.50p and (I2.5ip guarantee that a non-zero momentum 
is imparted to the particle after interacting with the plane gravitational wave. 
Hence gravitational waves carry energy. 
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2.4 Sources of gravitational waves 



2.4.1 Spinning rod 

A rod spinning about an axis perpendicular to its length is one of the first 
sources of gravitational radiation ever to have been considered [13]. 

Consider a steel beam of radius r = Imeter length / = 20 meter, density 
p = 7.Sgm/cm^, mass M = A.dxlO^gm and tensile strength t = 40, OOOlbft/in^. 
Let the beam rotates about its middle, so it rotates end over end with an an- 
gular velocity uj limited by the balance centrifugal force and tensile strength. 

a; = (-^) V2 = 28r adians/ sec (2.52) 
The internal power flow is 



Linternal = i^^'^'^)'^ = ^^^^'^^ 

2 X lO^erg/sec ^ 10'^^ Lq joule (2.53) 

where Lq = c^/G. The order of magnitude of the power radiated is only 
Law ~ {lO~^^yLo ~ 10~'^^erg/ sec . Evidently the construction of a laboratory 
generator of gravitational wave is unattractive, as this is a very small quan- 
tity which cannot be detected without new engineering or new ideas. There 
are however a great variety of astrophysical sources of gravitational waves. 
We list some of them and then discuss them lightly without going into actual 
calculation. 
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2.4.2 Astrophysical sources of gravitational waves 

Astrophysical sources of gravitational waves are the following : 

1. (a) Pulsar 

(b) Double star system; 

(c) Gravitational collapse of a few solar mass star; 

(d) Formation of a large black hole;. 

Now let me explain them. 

(a) Gravitational radiation from a pulsar 

Consider a highly dynamic astrophysical system. In particular take it to 
be a wildly rotating pulsar. If its mass is M and its size is R then by virial 
theorem its kinetic energy is ~ M"^ /R^ . The characteristic time scale for mass 
to move from one side of the system to the other is 

T ~ R/ (mean velocity) ~ R/iM/Rf^ = {R^/Mf^ 

(2.54) 

The internal power flow is 



HneUcenergy ^ (^)(^)i/2 ^ (^/^)5/2_ (2 55) 

T R R"^ 



The gravitational wave output is the square of this quantity or Lqw ~ {^)^Lq. 
Clearly the maximum power output occurs when the system is near its gravi- 
tational radius, and because nothing, not even gravitational waves can escape 
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from inside the gravitational radius. The maximum value of the output is 
~ Lq = 3.63 X IQi^^erg/ sec regardless of the nature of the system. 

(b) Double star system 

It has been estimated that at least one-fifth of all the stars are binary 
systems. We will go into the details of how it happens so frequently because 
that is a topic of astrophysics and hydrodynamics. 

Consider two stars of masses mi and m2 revolving in a circular orbit about 
their common centre of gravity. For their circular frequency of revolution, we 
have the standard formula: 

uj'^ — (mi+m2)/r^ (geometrical units of mass and time) (2.56) 

The calculated rate of loss of energy by radiation is [13]. 



- ^ = (^)(mim2/(mi + m2)) W (2.57) 
at o 

The following are important types of double star systems. 
Type name —{^)grav on Earth 

77 cos 5.2 X IQi^^erg/ sec 1.4 x 10~^^erg'/cm^ sec 

il3oo 3.6 x lO^^ergf/sec 6.7 x lO'^^ergi/cm^ sec 

Wuma 4.7 x lO'^^erg/sec 3.2 x 10~^^erg'/cm^ sec 

(c) Gravitational collapse of a few solar masses star 

Collapse to form neutron stars or black holes in the mass range 1 to lOM© 
{Mq = 2 X lO^^gm is a solar mass) will radiate waves in the frequency range 
1 to 10 kHz with an amplitude that depends on how much symmetry there 
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is in the collapse. These collapses, at least sometimes, result in supernova 
explosions. The rate at which supernova occurs is relatively well known, but 
the fraction of collapse events that produce strong enough gravitational waves 
is not well known. The characteristic period of the waves is proportional to the 
light-travel time around the collapsed object, the dominant frequency scale is 
1/M. For sufficiently large M, the source will produce low frequency waves 
detectable in space. 

(d) Formation of a giant black holes 

Many astrophysicist believe that the most plausible explanation for quasars 

and active galactic nuclei is that they contain massive ( 10^ — IO^Mq) black 

holes that accrete gas and stars to fuel their activity. There is growing evidence 

that even so called normal galaxies, like our own and Andromeda, contains 

black holes of modest size (10^ — IO^Mq) in their nuclei. It is not clear how 

such holes form, but if they form by the rapid collapse of a cluster of stars or 

of a single supermassive star, then with a modest degree of non-symmetry in 

collapse, they could produce amphtudes h fa 10~^^ to 10~^®meters in the low 

frequency range observable from space. If a detector has spectral noise density 
_ 1 

of lQ~'^^Hz ^then such events could have signal to noise ratio (;^) of as much 
as 1000. This strong signal would permit a detailed study of the event. If 
every galaxy has one such black hole formed in this way , then there could be 
one event per year in a galaxy. If no such events are seen, then either giant 
black holes don't exist or they form much more gradually or with too much 
spherical symmetry. 
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Chapter 3 

REVIEW OF THE 
PSEUDO-NEWTONIAN 
FORMALISM AND ITS 
EXTENSION 

In the first section a review of the Pseudo- Newtonian (ipN) formahsm is given. 
Section 2 provides a review of the extension of the ipN formahsm. In section 
3 the extended t/jN [eipN] formalism is used to develop a formula for the 
momentum imparted to test particles in arbitrary spacetime. Finally this 
formula is applied to plane and cylindrical gravitational waves, both of which 
give very reasonable results. 
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3.1 The ^A^ formalism 

The formalism [7] is based on the observation that, whereas the gravita- 
tional force is not detectable in a freely falling frame (FFF), that is so only at 
a point. It is detectable over a finite spatial extent as the tidal force. It could 
be measured by an accelerometer, as shown in fig.3.1. 

This accelerometer has a spring of length "/" which connects two masses. 
The spring ends in a needle which can move on the dial of the accelerometer to 
give a measure of the tension in the spring. Thus an observer in the FFF can 
observe the position of the spring by observing the moment of the needle on 
the dial. In Newtonian gravitation theory, the needle will show a zero position 
of the accelerometer in the absence of a central force. The force exerted by 
the source pulls the mass near it more than the mass further away. Thus the 
spring is stretched and the needle will move in the positive direction. 

If the spring is compressed the needle will show a negative deflection oth- 
erwise it will show a positive deflection. This would occur if both discs and 
the source consisted of like electric charges. Hence the negative deflection cor- 
responds to a repulsive source and the positive deflection corresponds to an 
attractive source. The strength of the source would be shown by the extent 
that the needle moved. 

Mathematically the tidal acceleration is given by 



41 



A^ = R^,^J''Fe (/x,z/,... = 0,1,2,3) 

(3.1) 

Here Ris the Riemann tensor, t a timelike killing vector, 1 is the spacelike 
separation vector representing accelerometer. Thus in geometrical terms tidal 
force is given by 



F'^ = -mR^pJ^Fr (3.2) 



where m is the mass of the test particle, the timelike vector tan- 
gent to the particles path and the separation vector, which provides the 
observation of the tidal force. In the FFF = /5q, P ~ Thus equation 
(13.21) becomes: 

= -mpR^opol'- (3.3) 

Regarding this as an eigenvalue equation, we get the maximum tidal force 
along the eigenvector of the matrix R^pQ. Since F is a purely spacehke vector 
in the free fall frame, the maximum tidal force will be a purely space like 
vector. Thus we get 



fWt:, (z,j = 1,2,3) (3.4) 



where " ; " stands for covariant derivative. Thus 7* is the relativistic analogue 
of the Newtonian gravitational force. Further this force is the gradient of a 
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Figure 3.1: An accelorometer to measure the force due to a gravitating source. 
The force is measured by the instrument in terms of the tension in the spring 
and shown by motion of the needle on the dial. 
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scalar quantity [7] 



(3.5) 



Where an appropriate expression for (f) is 



0= 2(^00-1)- 



(3.6) 



When expressed in this way a Lorentz factor has to be introduced by hand. 
7' and are called the tpN force and potential respectively. 



The quantity whose directional derivative along the accelerometer, placed 
along the principal direction, gives the extremised tidal force, which is zero 
in Minkowski space. Thus the e-^tN force, F'^, satisfies the equation [9] 



Where F*'* is the extremised tidal force. The fact that the zero component of 
the left side is zero does not guarantee that the zero component of F'* is zero. 

The space could be made look fiat in a small neighborhood of any point P 
in by a special choice of inertial coordinates. This locally inertial coordinate 
system has the metric of special relativity at some point P and so the Christof- 
fel symbols F^^ are zero at P. It is not, however, possible to make RfxvpK = at 
P by any choice of coordinates, unless the space is flat in the neighborhood of 



3.2 The ei)N force 



(3.7) 
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p. This is because R^pp-n is a tensor; if it vanishes in one coordinate system it 
will do so in any other because of the transformation laws. Here we will restrict 
ourselves to Riemann normal coordinates spatially, but not temporally. 
Equation (13.71) can be written in the space and time break up as 

r(FO + r°F^) = (3.8) 



F(F;. + r^j.F") = (3.9) 

A simultaneous solution of the above equations can be obtained by using Rie- 
mann normal coordinates for the spatial directions but not for the time coor- 
dinates. Thus the eipN force four vector is (see Appendix 2 for the proof): 

F° = m[(ln/l),o - + Tl^nj A]f\ (3.10) 

F^ = r,,f. (3.11) 
where A = (In y/—g)fi g = det{gij) For these, block diagonalized metrics 

J- 00 — 2 5' 5'00,0 , 00 — ~2fi^ 9oo,j 

Thus, writing the covariant form of the eipN force 

Fo = m[(ln A/),o - g%,fl/^A], (3.12) 
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F, = m(ln/),. (3.13) 

It is worth considering the significance of the zero component of the e%l)N 
force, which is its major difference from the ^l)N force. In special relativistic 
terms, which are relevant for discussing forces in a Minkowski space, the zero 
component of the four vector force corresponds to a proper rate of change of 
energy of the test particle. Further we know that, in general, an accelerated 
particle either radiates or absorbs energy according as ^ is less than or greater 
than zero. Thus Fq, should also correspond to energy-emission or absorption 
by the background spacetime. This point will be discussed further in the next 
chapter. 

3.3 General formula for the momentum im- 
parted to test particles in arbitrary space- 
time 

There has been a debate whether gravitational waves really exist [3,4]. To 
demonstrate the reality of gravitational waves Ehler and Kundt [4] considered 
a sphere of test particles in the path of plane fronted gravitational waves and 
showed that a constant momentum was imparted to the test particles. This 
was latter extended by Weber and Wheeler [3] for cyhndrical gravitational 
waves. The %l)N force and e%l)N force are based on an operational procedure 
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embodying the same principle [6]. The proper time integral of the force four 
vector will be the momentum four vector. Here it will be verified that this 
procedure gives the Ehler-Kundt result for plane fronted gravitational waves. 
When it is applied to cylindrical gravitational waves it is found that the result 
so obtained is physically reasonable and gives an exact expression for the mo- 
mentum imparted to test particles, corresponding to the approximation given 
by Weber and Wheeler. 

From Equations f l3.10p and (13. lip the momentum four vector, P^, is ob- 
tained as: 



Now we apply this formula to plane fronted gravitational waves and cylindrical 
gravitational waves. 

3.3.1 Plane-fronted gravitational waves 

The metric for Plane-fronted gravitational waves is 

ds"^ = dt^ - dx^ - L^{t, x){exp(2/3(t, x))dy'^ + exp(-2/3(t, x))dz^}. (3.15) 

Where L and (3 are arbitrary functions. 
The metric tensor is 




(3.14) 
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Its inverse is 



gab — 



1 



0-10 



-L^e^^ 



-L2e-2/3 



(3.16) 



9 



ah 



1 



0-10 



-L-2e-2/3 



V 



(3.17) 



-L-^e^^ 
To find the momentum four vector we need the force four vector. We calculate 



the force four vector. Here in this special case: 



g = dei{gij) = -L^(i,x), / 



(3.18) 



lny^ = 21nL 



(3.19) 



which implies that 



^/=(ln^),o=2^. 



(3.20) 



(ln(A/)),o = f - I 



(3.21) 
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8{(3 +L). 



(3.22) 



So from equations ^?M . ^?m . (13:221) and fl3:Tnll we get 



Fo = m{L + (3 L)/L. 



■ 2 



(3.23) 



From the vacuum Einstein equations we have (already discussed in chapter 2) 



.which implies that Fq = that is the zero component of the force four-vector 
is zero. Also Fj = because (In^/^) ^ = 0. 

Thus the momentum four vector becomes = constant. Hence there is a 
constant energy and momentum imparted to the test particles. The constant 
here determines the strength of the wave. This exactly coincides with the 
Ehler-Kundt method in that they demonstrate that the test particles acquires 
a constant momentum and hence a constant energy from a plane gravitational 
wave. 

3.3.2 Cylindrical gravitational waves 

Consider the cylindrically symmetric metric (12. ip . We first calculate the force 
four vector. Here in this special case: 



L + i3 L = 0, 



■ 2 



(3.24) 




(3.25) 
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In = 2(7 - ip) 



(3.26) 



A^2{r-t), f^e^-\ (3.27) 

(in(^/)),o = ^^^fy- - (r - 1). (3.28) 

ygy,,- = -4(7- - tf - (3.29) 

where a dot denotes differentiation with respect to t. 
The vacuum Einstein field equation gives [3]. 

ip{t, p) = AJq{up) cos{ut) + BNo{ujp) sm{ujt), (3.30) 



where Jq{ujp) and Nq{ujp) are the zero order Basel and Neuman functions 
respectively. A and B are arbitrary constants corresponding to the strength 
of gravitational wave. 

7(i, p) = ]^up{{A^Jo{up)J'^{up) - B''No{x)N;,{x)) cos{2ut) 
-AB[{Jo{ujp)N'^{ujp) + J'^{x)No{x))sm{2ujt) - 
-2{Jo{ujp)N'^{up) - J'^{x)N^{x))u:t]} (3.31) 

Where prime now refers to differentiation with respect to ujp, uj being the 
angular frequency. 
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Using equations (13.281) , (13.291) , (I3.30p and (13.311) we get the zero component 
of the force four vector. 



Fo = -muj{[AJo cos(u;t) + BNq sin(cjt)] - 2pcj[(AVo Jq - B^NqNq) cos{2ujt) 

, 2\AJo sm( cut) - BNq cos(iut)]'^ 

—BNq cos(co't) — 2ujp[AJo sin(ci;t) — BNq cos(ci;t)] 

[AJq cos(c<;t) + BNq sin(cjt)]}. (3.32) 
As 7 and ip are functions of t and p only so F2 and F3 are zero. 

Fi = m{\nf)^i = -muj{[AjQCos{iJt) - BN^siniut) -^[{A^JqJq- B^NqN;^) 
+uop{A^JqJq - B'^NqNq)'] cos{2ut) - ]^AB[2{JqNq + J'qNq)'] sm{2uot) 
-]^AB[A{JqN'^ - J'^Nq) + 2ujp{JqN'^ - J^iVo)Vt}. (3.33) 

The corresponding Pq and Pi are 
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Pq = — m[ln I A Jq sin(co't) — BNQCOs{ujt) | + In | 1 — 2ujp[AJ'QCOs{ujt) 



_i (1 + 2AojpJ'q) tan(|cjt) - 2BupN'^ AB{NqJ'q - J^N'^ 

^""'^ ' - 4a;2p2(^2 JV2 + 52Ar^2) I +p(A2 + ^^^^^2)^ 

+/i(o;p).(3.34) 

Pi = -m{[AJo sm(c<;t) - PA/'q cos^)] - ^[(AVq Jq - B^NoN^) 
+ujp{A^JoJo - B^NoK)'] sm{2ujt) - ^^^[(JoiV^ + J^iVo) 
cup(JoiVo + ^o^o)'] cos(2cut) - A5a;'t2[(JoiVo + 4^o) 
-a;p(JoiVo - W)']} + /2M- (3.35) 



Where /i and /2 are arbitrary constants of integration. Weber and Wheeler 
exclude solution that contain irregular Bessel function, Nq^up), as not well 
defined at the origin. Taking the Weber- Wheeler solutions equations fl3.34p 
and (I3.35P reduces to 



Pq = -m[ln I y4Josin(tut) I +(1 + y4Jo/cjpJo)ln I 1 - 2cjpAJoCos(cj03|36) 




52 



Pi = -m{AJ^ sin(cjt) - -A'^JqJq + pcu( Jq Jq)' sm(2cjt) + /2(c^p)}. (3.37) 

We see that the quantity Pi given by equation (13.371) can be made zero for the 
large and small p limits by choosing /2 equal to zero. This is physically rea- 
sonable expression for the momentum imparted to test particles by cylindrical 
gravitational waves. The quantity Pq given by equation (I3.36P remains finite 
for small p and can also be made finite for large p by choosing /i = — In(Jo). 
However, there is a singularity at ut = nn. This problem does not arise in 
the general expression given by equation (I3.34p . However, in that case there 
appears a term linear in time which creates interpretational problems. Also 
Pq and Pi becomes singular at p = if P 7^ 0. 
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Chapter 4 



SPIN IMPARTED TO TEST 



PARTICLES BY 



GRAVITATIONAL WAVES 



We can write the e'il)N force four vector (discussed in section 3.2) as [7] 



i^o = -C/,o , Fi^-Vi. (4.1) 



where 



U = m 



(4.2) 



y = -mln(/), (4.3) 

where B is a constant with units of time inverse so as to make ^ dimensionless. 
Here V is the eipN potential but there is no good interpretation of U. Thus 
there is a problem of interpretation of Fq. The Fi was reinterpreted [8] as the 
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rate of change of the momentum imparted to test particles by the gravitational 
field, i.e. Fi = ^ (where r is the proper time). In this interpretation it would 
be natural to identify Fq as The problem now is to interpret Pq since one 
would normally take Pq = E = ^ + P^P^gij where m is the mass of the 
test particle . Hence / Fodt cannot be this Pq. Sharif's suggestion [11] for 
the interpretation of Pq, that it gives the spin angular momentum imparted to 
test rods, is given in the first section of this chapter, but in the same section 
it turns out to be inconsistent. To find an alternative check on its validity the 
geodesic [4] analysis for the angular momentum imparted to test particles by 
gravitational waves is undertaken in section 2. This formalism is applied to 
various cases in section 3. 

4.1 Spin angular momentum imparted by grav- 
itational waves 

Sharif [11] considers a test rod of length A in the path of a gravitational wave 
whose preferred direction is given by /* in the preferred reference frame. He 
argues that the rod will acquire maximum angular momentum from the wave 
if it lies in the plane given by epjkil'^, where e^,^p,r is a totally skew symmetric 
fourth rank tensor. Thus the spin vector will be given by 

= le'""'ejkil'P.. (4.4) 
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For i — 1, 



S 



1 







= Pol'. 



(4.5) 



Similarly 



Pof and = Pol^. 



(4.6) 



So in the preferred direction the spin vector would be proportional to T such 
that 



The angular momentum imparted would be the magnitude of the spin vector 
. Thus the maximum angular momentum imparted to a test rod when it lies 
in the plane perpendicular to the preferred direction is: 



Hence the physical significance of the zero component of the momentum four 
vector would be that it provides an expression for the spin imparted to a test 
rod in an arbitrary spacetime. 

This formula was applied to plane and cylindrical gravitational waves to 
give the following results. 



Pol'. 



(4.7) 




(4.8) 
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1. Plane gravitational waves. 

Using metric (12.251) in equation (14. 8 p we get (for detailed calculations see 
section 3.2) 

Pq = constant. (4.9) 

and thus the spin would also be constant. 

2. Cylindrical gravitational waves. 

Following the same procedure for the metric (2.1) we get (for detailed 
calculations see section 3.2) 



S = — mA[(l + AJq/ojpJq) In | 1 — 2uj pAjQCOsiut) \ +constant. (4.10) 

Notice that there can be no spin angular momentum imparted to test par- 
ticles in a perfectly homogeneous and isotropic cosmological model [1]; its high 
degree of symmetry — in particular, spherical symmetry is incompatible with 
spin being imparted to test particles. However when we use Sharif's formula 
for cosmological models, it gives exactly this error. 

We give examples which had already partly been constructed by M. Sharif 
[14]. 

3. The Friedman model: 

Consider the Friedman model, which is isotropic and homogeneous. 



ds' = dt' - a\t)[dx' + fMdn% 
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(4.11) 



where 

k = l = sin(x), 

k = foix) = X, 

k = -l = sinh(x), 

X being the hyperspherical angle and a{t) the scale parameter given by: 
a{t) = ao(l — cosri)/2, t = ao{ri — sin//), < f] < 2n, k = 1 



a, 



= ao(cosh?7 — l)/2, t = ao^smhr] — ri)/2, < r/ < oo, k = —1; 

(4.12) 

Using the metric (14. lip the eipN force for the Friedmann models will be 

Fo = -m-, (4.13) 
a- 

and 

Pq = J Fodt = -mln(a ). (4.14) 
As S = mXj Fodt. We get the spin for Friedmann models 

S = -m\\n{a-). (4.15) 

a. Closed Friedmann model 

Using the metric (14. lip with = 1 in equation (14.140 and (I4.15P we get 

Pq = m[ln a/1 — COST] cos?7 H cos^ 77 + In y/2 H ], (4.16) 

8 16 16 



and thus 



3 ]^ 5 

S = mX [In \/l — cos rj — - cos V + Tc "^^^^ rj + \n \/2 + — ] . (4-17) 

8 16 16 
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b. Flat Friedmann model. 

Again Using the metric (14.111) with /c = in equation (14.141) and (I4.15P we 

get 

Po = mln(|). (4.18) 
And thus the spin for a flat Friedmann model is 



^ = mAln(^) . (4.19) 



.2 

c. Open Friedmann model. 

Finally we obtain the spin for the open Friedmann model by the same 
procedure as the following 

Po = -mln| I- (4-20) 



cosh T] — 1 

sinhr/ 
cosh Tj — 1 



S = -mXln\ |. (4.21) 



Equations (I4.17P , (I4.19P and (I4.2ip tells us that a non-zero spin is imparted 
by flat, open, and closed Friedmann models. Since no spin can be imparted 
thus Sharif's interpretation of Pq cannot be correct. 

4. The Kasner model 

Consider the Kasner model for a homogeneous anisotropic universe (near 
the cosmological singularity) 

Here pt are numbers such that pi + P2 + Ps = Pi + P2 ~^ P3 = ^■ 

59 



The metric tensor is 



Its inverse is 



( \ 

10 



-t^pi 







y -t^ps J 



1 











v 



(4.23) 



(4.24) 





To find the momentum four vector we need the force four vector. We 



therefore calculate it. Now 



g = det{gij) = -t^ / 



In ^/—g = In t, 



(4.25) 
(4.26) 



which implies that 



A/=(lnv^),o 



(ln(A/)),o 



I] 

9,o9ij,o 



So from equations (14.281) and fl4.29p we get 



(4.27) 

(4.28) 
(4.29) 



(4.30) 
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Hence from equation (14. 8p we get 



5* = constant. 



(4.31) 



By physical consideration we could set 5 = 0. 

5. The De Sitter universe (usual coordinates) 

Consider the metric 

ds^ = (1 - r^/D^)de - (1 - r^/D^)-^dr^ - r^dO^ + sin^ Odip''), (4.32) 

The metric tensor is 

^ (l-rVD2) 













[l-r^D^y^ 
—r 



\ 














-r^ 
-r^ sin^ 9 



Its inverse is 



\ 










-4 













(4.33) 



(4.34) 



Again proceeding on the same lines for the momentum four vector, we have 



g = det{gij) = -r sin 6, f 



In \/—g = ln(r^ sin ( 



g^ v^(l -rVD2) 



(4.35) 
(4.36) 



which implies that 
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A/ = (lnv/^),o = 0, 



(ln(A/)),o = 0, 



9fl9i3,o = 0. 
So from equations (14.381) . (I4.39p we get 



(4.37) 
(4.38) 
(4.39) 



Fn = 0. 



Hence from equation (14. 8 p we get 



S = constant. 



(4.40) 



(4.41) 



Clearly, we would need to take 5" = here. 

6. The Lemaitre form of the De Sitter universe 

Consider the empty space solution of the Einstein field equations with 
cosmological constant, 

ds^ = de - a2e2(^/3)^/'* [dx^ + x^de^ + sin^ Od^^] . (4.42) 



The metric tensor is 
1 



9fj.u 



_a2e2(A/3)V2i 











^2g2(A/3)V2j^2 gi^2^ 



(4.43) 
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Its inverse is 



\ 









\2 JVsVAt 













(4.44) 



To find the momentum four vector we need the force four vector. We 



therefore calculate it. Now 



/goo 



(4.45) 



In = 3 In ao+2 In x+v^VAt+^ In (- cos ^ + 1)+^ In (cos ^ + 1) , (4.46) 
which implies that 



A/=(lnv^),o = v^VA 



(4.47) 



(ln(A/)),o = 0, (4.48) 



9%,,o = -4A. (4.49) 



So from equations fICTD . flOOj) and fl332D we get 



A 



1/2 



Fo = m[-] . (4.50) 



Hence from equation (14. Sp we get 



A 



1/2 



S = mX ( ) ^ + constant. (4-51) 
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This does not seem reasonable. Since the other form of the metric gives a dif- 
ferent result, it is clear that the interpretation is not even internally consistent. 

4.2 The geodesic analysis for angular momen- 
tum imparted to test particles by gravita- 
tional waves 

In section 4.1 we have concluded that S is not the spin angular momentum 
imparted to test particles by gravitational waves. Then the question arises, if 
S is not the spin then what is the spin angular momentum. So in this section 
we use the geodesic analysis [4] to find the angular momentum imparted to 
test particles by gravitational waves. This formula is further applied to various 
cases. 

Consider a time like congruence of the world lines (not necessarily geodetic) 
with tangent vector rt". Decompose Ua-b by means of the operator hab projecting 
in to the infinitesimal 3-space orthogonal to [4]. 

Ua-b = -^^ab + CTab + \oKb " U'^Ub. (4.52) 

where 

- UJab = U[a,b] + Ul^Ub] , = Ua-fiU^ ■ (4.53) 

= u'^a > f>'ab = gab + Ma^ft (4.54) 

CTab = U^a;b) + ^^(a^^b) " ^dhab « = 0). (4.55) 
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For an observer along one of the world lines and using Fermi propagated axes, 
ujab describe velocity of rotation, aab shear and 9 describe expansion of the 
neighboring free particles. Since the ■?/'A^— formalism uses the fermi-walker 
frame, it could be expected that the results of this analysis should be consistent 
with it. For our purpose only Uat is needed. Choose the coordinates so that 
the tangent vector is = /5q. In the case goo = 1, = 6q. If goo 7^ 1, = 
-^^^o Thus, from equation fl4.53p we have 

- Uab = {Ua,b - Ub,a)/'^ + {UaUb - %Ua)/2. (4.56) 

The first and second component on the right hand side of the above equation 
vanishes. Thus we have 

- uJab = {u-gUb - u-f;aa)/2. (4.57) 
Further using equation (14.531) we have 



2uj. 



ab 




Taking a to be i and h to be zero, we finally obtain the components of the spin 
vector: 

^,o = -^<! " )■ (^ = 1,2,3) (4.59) 
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This simple formula appears because u"' — ^Qonly. This gives the angular 
momentum imparted to test particles by gravitational waves. We now apply 
this formula to different types of gravitational waves. 

4.3 Applications 

4.3.1 Plane gravitational waves 

For a metric gij_i, construct the new metric 



(4.60) 



where 77^1^ is the Minkowski spacetime and u{s — sq) is the step function defined 



as: 



u{s — So) = when s < so 

u{s — So) — 1 when s > Sq . 

.According to this definition the metric tensor becomes: 

^ [(e2^(") - 1)m + 1] 

(1 - e^^^°''^)u - 1 

(-Q;2e2/3(") + l)u-l 









_a2g-2/3(a) ^ ;L)« - 1 j 

(4.61) 



where a — t — x. The inverse of this metric is 
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(e2f2{")-l)M+l 










v 



Now 





1 



^l+(-e2"(")+l)n 





25'°°fl'oo,i5 









1 

-l+(-a2e2/3(«)+i)„ 







2 



25'°°fl'oo,2, 



-l+{-a2e-2'3(")+l)u / 

(4.62) 



3 



25'°°5'oo,3- 



(4.63) 



If we find gm,ii (700,2 and (yfoo.s then we are done. Since qqq = (e — 1)^ + 1. 



(4.64) 



ds 



2fi'(«)e2^(")M - (e2^(°) - l)^5(s - sq). 

ox 



Dividing equation (12.251) by ds'^ and keeping t,y and z fixed, we get: 



9s 
dx 



(4.65) 



Consequently we have 



5'oo, 



1 = -21]'(a)e2^(")M(s - So) + e^^^^e^^^") - l)6{s - Sq). (4.66) 



Further 



^700,2 



3^^(")-l)|^^(^-^o). 



(4.67) 
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Again dividing equation fl2.25p by ds^ and keeping t, x and z fixed, we get 

|^ = ae^(°). (4.68) 
dy 

Putting this value of|| in to equation (I4.67p we get 

^700,2 = -e^(")(e2^(") - - sq). (4.69) 

Also, 

^7oo,3 = -«(e^^(")-l)|^5(s-so). (4.70) 

As in the previous cases, dividing equation fl2.25p by ds^ and keeping t, x and 
y fixed, we get 

|^ = ae-^("). (4.71) 
oz 

Putting d^ in equation (I4.70p we get: 

^?oo,3 = -ae-^(°)(e'^^"^ - l)5(s - So). (4.72) 
Using equations (14.660 . (I4.69P and fl4.72p in equation (I4.63P we get: 





1 


2 



1 



2(e^^W - l)n + 1 [-^^^^^'''Ms - So) + e-(")(e^-(") - 1)6(8 

(4.73) 

(4.74) 
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(4.75) 

Now using equations (14. 73 p . f l4.74p and (14.751) in equation (14.591) we get com- 
ponents of the spin vector 




^10 = - 4(,2n(.)_ 1)^ + 1 [-2Q'ia)e^^^Ms - So) + e-(")(e^-(") - 1)6(8 - .o)] 

(4.76) 

(4.77) 

-30 = - 4(,2.(.)li), + i [«e-^(")(e-(") - ms - .„)] • 

(4.78) 

In equation (I4.76P the first part in the brackets gives the spin imparted to test 
particles by the plane gravitational waves while the second part gives the spin 
of the wave itself. 

4.3.2 Cylindrical gravitational waves 

Now use the same construction for the metric (12. ip to get 



^700 = (e2(^-'^) - l)w(s - so) + 1, 

922 = -p'[{e-'^~lHs-So) + l], 
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5^33 



{e'^ -l)u{s-so) + l. 



(4.79) 



The inverse of this metric is 



( on ^ 






^33 



(4.80) 



where 



00 



1 



(e2(7-i/')-l)u(s-so)+l ' 

11 _ 1 

_(e2{7-'^)_l)„(s-so)-l' 

22 _ 1 

-p2[(e-2^-l)n(s--so)+l] 
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1 



-(e2*-l)tj(s-so)-l ■ 

In this case 



^00,1 = 2e2(^-^)(y - ^')u{s - So) + (e2(^-'^) - l)^^6{s - s^y (4.81) 



Dividing equation (12.11) by ds"^ and keeping t, (p, z fixed, we get 



1 _ .2(7-^) {^\ 



(4.82) 



which imphes that 



dp 



(4.83) 



Now putting 1^ in equation (14.811) we get 
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^700,1 = +2e^(^-'^) (7' - u{s- so) + e^"'^ {e'^'-^^ -1)5(3- so) , (4.84) 



5^00 



^ = _[(e2(7-^)_ 1)^ + 1]^ 



3'^^-^^ -1)^5 (3 -So), (4.85) 



900,3 = - {e'^'-^^ - 1) {s - So) . 



Dividing equation (12.11) by ds'^ and keeping t,p, z fixed, we get 



(4.86) 



l = e2V(^)', (4.87) 



which imphes that 



= Pe-^- (4.88) 

dy 

Now dividing equation (12. ip by (island keeping t, p, y fixed, we get 

(4.89) 

l^e*. (4.90) 
Putting equation (14.881) in to equation (14.851) we get, 

9oo,2 = -{e'^^-^^-l)pe-^5is-So). 

(4.91) 
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Putting equation (14.901) in to equation fl4.86l) we get, 
^700,3 = - {e"^'-^^ - 1) eH {s - s,) . 

(4.92) 

Using equations (14.841) . (I4.9ip and (I4.92p in equation (14.631) we get: 





1 



[2e'(^~'^\j' - ij')u + e^-'^(e2(^-^) - 1)5(3 - so)] 



2(e2{7-^) _ 1)^ 

(4.93) 





2 


3 



(4.94) 



(4.95) 

Now using equations (I4.93p . (I4.94p and (14.950 in equation (I4.59P we get the 
components of the spin vector 

'-lo = ,(^.(,_,)'_,)^^, [2e^(-^)(y - ^')u + e-ne^^-'^) - 1)5(. - .o)] , 

(4.96) 

(4.97) 

(4.98) 
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4.3.3 The Friedmann model 



Consider the Friedmann model. The metric for this model is given by equation 
(14. lip and the metric tensor is 



\ 



(4.99) 



10 

-a2(t) 

-a\t)fi{x) 

-a\t)f,ix)smH 

According to equation (14.591) the only thing we need is 5^00,1, 5^00,2, and (700,3- 
As there are no wave fronts. Therefore we have Qoq = 1. Here obviously all it's 
derivatives will vanish. Hence 



v 



f 1 


hi 


\ 1 


hi 


I 1 


[o.J 











0. 



(4.100) 



Thus all components of spin vector are zero. i.e. 

U^io = iO20 = ^30 = 0. 



(4.101) 



As we have stated earlier, there is no spin in an isotropic and homogeneous 
universe model, this analysis gives zero spin for this model as required. 
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4.3.4 The Kasner model 



Consider the Kasner model. The metric for the Kasner model is given by 
equation (14.22p and the metric tensor is 



\ 



(4.102) 



10 

-t^Pi 


-t2p3 

According to equation (14.591) the only thing we need is (700,1 , ^700,2, and (700,3- As 
there is no "wave-front" so the original metric stands. From (I4.102p we have 
(7oo = 1- Here obviously all it's derivatives will vanish. Hence 



V 



f 1 


hi 


\ 1 


hi 


\ 1 


[o.J 











0. 



(4.103) 



Thus all components of spin vector are zero. i.e. 

LOlQ = LO20 = 1^30 = 0. 



(4.104) 



4.3.5 The De Sitter universe (usual coordinates) 



The metric for this model is given by equation (14.321) and the metric tensor is 



;i - r^/D^) 









\l-r^/D^)-^ 

-r 








-r^ 
-r^ sin^ 9 



(4.105) 
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As goo 7^ 1 so ujio calculated in section 4.2 cannot work. Let u°- = /5q. Here 
by definition 

u'^u'^gab = -1 = (m°)2(7oo or m° = -L (4.106) 

Vfi'oo 

consequently we have 

= -^6^. (4.107) 

\/goo 



Now we only need to calculate u^a^b] and MLWh] 



«a,6 = (v^^a),i, 

= (lny^),i5>i. (4.108) 



Hence we have 



«M = (lnv/^),i5pX] (4.109) 



Similarly we get 



1 







u 



[a^l>] = -7^{Sya-Syb)V9^0{ } (4.110) 

1 

From equation ( 14.109^ and (14.110p the sum U[a,b]+u'^^Ub] for this metric vanishes 
and thus equation (14.531) implies that ujab vanishes. 
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4.3.6 The Lemaitre form of the De Sitter universe 



The metric for this model is given by ( ]4.42p . The metric tensor is 



9^lv 



\ 



_a2e2(A/3)^/^t 



















-a2e2(^/3)^/'*x' sin^ 9 



J 

(4.111) 



By the definition (I4.60p of the metric tensor qqq = I.Here obviously all it's 
derivatives will vanish. Hence 



f 1 


hi 


\ 1 


hi 


\ » 1 


1-J 











0. 



(4.112) 



Thus all components of spin vector are zero. i.e. 

tUiO = UJ20 = ^30 = 0. 



(4.113) 



Thus we have got the same result for both the forms of the De Sitter universe 
which proves that the formalism is consistent. 

4.3.7 The G odel universe model 

Consider the Godel (spinning) universe model 

ds"^ = dt^ - dx^ - dy"^ + 2m{x)dtdz - l{x)dz^ , (4.114) 
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where m(x) = Ae^^ and l{x) = ^^(tj — l)e^'^^ . The metric tensor i 



IS 



1 m(x) 
0-10 

0-10 

m{x) 



\ 



(4.115) 



Its inverse is 



l{x) 



l[x)+m'^ (x) 




mix) 





-1 

-1 





•m(x) 



\ 



l{x)+m? {x) 






1 

' l{x)+m?(x) J 



(4.116) 



y l{x)+rn?(x) 

According to equation fl4.59p the only thing we need are the following three 





Christoffel symbols 



1 


1 



f 1 




\ 1 






> and < 




> . Now 


1-J 









m[x) 



2 l{x) + m?{x 

1 arn?{x) 

2 l{x) + m?{x) 



-(Ae""),i 



f 1 


hi 


\ 1 


[o.J 







0. 



(4.117) 



(4.118) 



Thus all components of the spin vector are 

1 am?{x) 



1^20 



4 l{x) + m?{x) 
= 0. 



(4.119) 
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Chapter 5 

SUMMARY AND 
CONCLUSION 

Linearized General Relativity predicts gravitational waves. These waves are 
analogous to electromagnetic waves. There may be different types of gravi- 
tational waves i.e. plane and cylindrical gravitational waves etc. This back- 
ground was reviewed in chapter 1. Work on plane and cylindrical gravitational 
waves [3] has been helpful in understanding them further. The question of the 
reality of gravitational waves was discussed in chapter 2. There we used Weber- 
Wheeler [3] method for a particle in the path of a plane gravitational wave and 
obtained the standard constant momentum imparted to the particle. Some as- 
trophysical sources of gravitational waves were also discussed. In chapter 3 
the work of Qadir and Sharif [8] in which a general formula is developed for 
the momentum imparted to test particles by gravitational waves in arbitrary 
spacetime, was reviewed. In this paper the problem of the identification of 
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their zero component of the momentum four vector was mentioned. Sharif 
[11] gave a suggestion that it could be interpreted as the spin imparted to a 
test rod in an arbitrary spacetime. This suggestion was reviewed in chapter 4. 
Further analysis by Sharif [14] was shown here to provide a counterexample for 
this suggestion i.e. it gives a non zero spin for an isotropic and homogeneous 
universe model. Further, when the example of the De Sitter universe was con- 
sidered, the interpretation gave different results for the static and Lemaitre 
forms. This proved that the interpretation is not even internally consistent. 
As such Pq cannot be interpreted as the spin imparted to test rods. A geodesic 
analysis [4] was used in section 4.2 to evaluate the spin imparted to test parti- 
cles in various cases. In the cases of plane and cylindrical gravitational waves 
we got very reasonable results. As required , the Friedmann, the Kasner and 
De Sitter models do not impart spin to test particles according to this analysis. 
Further we got the same results for both forms of the De Sitter universe, which 
confirms the validity of our analysis. Finally, the Godel universe model gave a 
non-zero (non-constant) angular momentum, precisely as it should. 

If Pq does not provide an expression for the spin imparted to test rods 
in an arbitrary spacetime then what is its correct interpretation?. Consider 
E = ^ P'^P^Qij -|- and Pq = J Fgdi. Define the difference between the two 
as AE — P^ — E for some of the spacetimes considered earher. By definition 
the usual energy E associated with a momentum is given by 

P'^P'^g^, = = {E^) goo + P'P'giy (5.1) 
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Hence 

E=^{m'-P^P^g,,f\ (5.2) 
Now, we want to compare this with 

po = J Fodt. (5.3) 

The energy difference, AE, is then 

= g'' [ Fodt - -L (m^ - P^P^g,^) . (5.4) 
J \'9oo 

When Fq is zero we get no further understanding from AE. As such we will 
not compute these cases. Again for Godel universe we have not calculated 
Fq as the metric is not in the block diagonalized form {g^i = 0) and a gauge 
transformation would need to be made before it could be applied. Therefore 
we only work out AE for the cases (a) cylindrical gravitational waves; (6) the 
Friedmann models and (c) the Lemaitre form of the De Sitter universe. Finally 
we will discuss the consequences of these results, 
(a) Cylindrical gravitational waves 

Consider the metric given by equation (12. ip and use the non-zero compo- 
nents of momentum four vector given by equations fl3.36p and fl3.37p to get 



^ n 2 \ 1/2 

AjQsm{cut) - -A^JqJq + pLU { Jo JqY sm{2ujt) 



+ 1 
(5.5) 

Consequently we have 



A^ = -me-(^~'^)[(l-e-2(^-^H^^oSm(cc't)-^AVoJo+pcj(JoJo)'sin(2u;t))2)i/2 

80 



+ e-(7-^)(ln I AJosm{ujt) \ +(1 + AJo/cupJo) In | l-2ujpAJ'QCos{ujt] |)]. 

(5.6) 



(b)The Priedmann models 

Since P^,P^ and are zero for the Friedmann models, from equations 
(14. lip . (14.171) . (14.181) and (14.201) we get the energy difference for the three 

models: 



SI 5 

AEk=i = m[ln \/l — cos rj — - cosr] + — cos^ rj + \n \/2 + —]; (5.7) 

AE,=o = mln(|); (5.8) 

A , coshr? — 1 , ^, 

AEk=-i = mln ^ . (5.9) 

sinh r] 



These are plotted in Fig. (5.1) , (5.2) and (5.3) and discussed shortly, 
(c) The Lemaitre form of the De Sitter universe 

Again P^, P^ and P^ are zero for the metric given by equation (I4.42p . Thus 
using equation (I4.5ip we get 

AE = m^/X/3t. 

(5.10) 

as shown in Fig (5.4) . 

The three expressions for the Friedmann models have the same asymptotic 
behavior for sufficiently small values of rj, namely AE ~ mlnr]. For k = 

2 

±1 there is a correction term ~ ='=12- have inserted a constant term for 



A; = +1 so that the expressions for the three should match up to the zero 
order terms. At the phase of maximum expansion of the closed model we get 
/S.E — m(ln2 + |). We could equally well, have set /S.E = at the phase of 
maximum expansion and had a difference for it from the other two cases for 
small 77 (in the constant term) . Note that for all the three models /^E diverges 
as 77 — > and it also diverges as 77 — > 27r for the closed model. For k — 
we have chosen to display the constant term so that /S.E = at 77 = tt. The 
proposal does not seem inconsistent, but still needs further discussion. 
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1 2 3 




Figure 5.1: Fig 5.1: AE for the closed friedmann model. This becomes infi- 
nite at the bang and the crunch and becomes zero at the phase of maximum 



expansion. 
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Figure 5.2: Fig 5.2: AE for the flat Friedmann model, again m is taken to be 
unity. AE diverges at r/ = 0. 



1 2 3 4 5 




Figure 5.3: Fig 5.3: AE for the open Friedmann model which diverges at = 
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OO 1 2 3j 4 5 6 



Figure 5.4: Fig 5.4: AE' for the Lemaitre form of the De Sitter universe. Here 
A is taken to be unity. 

APPENDIX 1 

Consider the metric [1] (with signature (+, — , — , — )) 



ds^ = L?{u){dx'^ — dy^) — dudv 



(5.11) 



where u = t ~ z andf = t + z which is always flat. 



The metric tensor is. 



9ah 







-L'^(u) 



-L'^iu) 







-1 



(5.12) 



Its inverse is. 
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^ I Q 



-L-'^{u) 

-1 
The non vanishing Christoffel symbols are: 



V 



(5.13) 



f ^1 


hi 


\ ^ 1 


hi 


\ ^ \ 


1-J 











L'{u) 
L{u) 



(5.14) 







[ 3 1 






hi 




> = L{u)L'{u) 











(5.15) 





l-l 




[a. J 







L{u) 



(5.16) 



1 

2 2 



3 
2 2 



= L{u)L'{u) . 



(5.17) 



[^1 


hi 


\ ^ 1 


hi 




[3. J 











L{u) 



(5.18) 



The primes refers to differentiation with respect to u. The nonvanishing com- 
ponents of the Ricci tensor are Rqq, Rq^, i?3oand Rss. Using Christoffel sym- 
bols listed above we obtain the following Ricci tensor components. 
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-2L"(u) 
L{u) 



R30 = i?03 = '^^rv ■ (5-20) 

L{u) 

Prom equations (Al.9) and (>11.10) the required result follows. 
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APPENDIX 2 

We shall prove here that 



F' = m[(ln A),o - r°o + KTyA]f, (5.21) 



and 



F' = (5.22) 
are the solutions of the following equations. 



f (F° + r°.F^) = 0, (5.23) 



P(F;. + r^^.F°) = (5.24) 

where 

F*' = m(no,, - r^,,o + K^Ko - nk^o)fl'- (5.25) 
For the verification of equation {A2.3) we note that 

A = f = — . (5.26) 

900 

Thus 

(In A),o = 9%l9'' + 9r3flo/9ij,o. (5.27) 

Here we will make use of the Riemann normal coordinates (RNCs) for spatial 
directions. We have 

9oi = (5.28) 

9tiy,o 7^ 0, g^y., = = g^yfii. (5.29) 



Using this approximation we can write 

(ln^),o, = = r^o,, = n,nM- (5-30) 

Now 

+m[{\nA),o - C + Tl^TyA]{f),. (5.31) 

Using RNCs the first, second, third and the last term on the right hand side 
vanishes. Therefore 

= (5.32) 

Thus equation (A2.3) is satisfied. Equation (^42. 4) can directly be obtained 
just by replacing the values of F° and from equations (^2.1) and (^42. 2). 

Hence F° and given by equations {A2.1) and (A2.2) are the solutions 
of equations (^2.3) and {A2A). 
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